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We construct a simple analytic model for wall-bounded turbulence, containing only four adjustable 
parameters. Two of these parameters characterize the viscous dissipation of the components of the 
Reynolds stress-tensor and other two parameters characterize their nonlinear relaxation. The model 
offers an analytic description of the profiles of the mean velocity and the correlation functions of 
velocity fluctuations in the entire boundary region, from the viscous sub-layer, through the buffer 
layer and further into the log-layer. As a first approximation, we employ the traditional return-to- 
isotropy hypothesis, which yields a very simple distribution of the turbulent kinetic energy between 
the velocity components in the log-layer: the streamwise component contains a half of the total en- 
ergy whereas the wall-normal and the cross-stream components contain a quarter each. In addition, 
the model predicts a very simple relation between the von-Karman slope re and the turbulent veloc- 
ity in the log-law region v + (in wall units): v + — 6re. These predictions are in excellent agreement 
with DNS data and with recent laboratory experiments. 

PACS numbers: 43.37.+q,47.32.Cc, 67.40.Vs, 67.57.Fg 



Introduction 



The tremendous amount of work devoted to under- 
standing the apparent experimental deviations from the 
classical phenomenology of homogeneous and isotropic 
turbulence [1, 2] tends to obscure the fact that in many 
respects this phenomenology is almost right on the mark. 
Starting with the basic ideas of Richardson and Kol- 
mogorov, and continuing with a large number of inge- 
nious closures, one can offer a reasonable set of predic- 
tions regarding the statistical properties of the highly 
complex phenomenon of homogenous and isotropic tur- 
bulence. Thus one predicts the range of scales for 
which viscous effects are negligible (the so called "in- 
ertial range" of turbulence), the cross over scale below 
which dissipative effects are crucial (also known as the 
"Kolmogorov scale"), the exact form of the third order 
structure function S^(i2) (third moment of the longitu- 
dinal velocity difference across a scale R), including nu- 
merical pre-f actors, and an approximate form of struc- 
ture function of other orders S n (R) (predicted to scale 
like i?™/ 3 but showing deviations in the scaling expo- 
nents which grow with the order, giving rise to much 
of the theoretical work alluded to above). In particu- 
lar much effort was devoted to calculate the so called 
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"Kolmogorov constant" C2 which is the pre-factor of the 
second order structure function, with closure approxi- 
mations (see, e.g. Refs. [3, 4]) coming reasonably close 
to its experimental estimate. Notwithstanding the de- 
viations from the classical phenomenology, one can state 
that it provides a reasonable first order estimate on many 
non-trivial aspects of homogeneous and isotropic turbu- 
lence. In contrast, the phenomenological theory of wall 
bounded turbulence is less advanced. In reality most tur- 
bulent flows are bounded by one or more solid surfaces, 
making wall bounded turbulence a problem of paramount 
importance. Evidently, a huge amount of literature had 
dealt with problem, with much ingenuity and consider- 
able success [5] . In particular one refers to von-Karman's 
log-law of the wall which describes the profile of the mean 
velocity as a function of the distance from the wall. It 
appears however that the literature lacks an analytically 
tractable model of wall bounded flows whose predictions 
can be trusted at a level comparable to the phenomeno- 
logical theory of homogeneous turbulence. 

In this paper we attempt to reduce this gap and offer 
as simple as possible but still realistic model (a "mini- 
mal model") for the viscous-to-turbulent flow in the en- 
tire region from the very surface through the logarith- 
mic boundary layer (hereafter, log-layer) up to the upper 
boundary of outer turbulent region. Our final goal is to 
create clear physical grounds for improved description of 
the flow-surface interactions in numerical fluid-mechanics 
models (both engineering and geophysical) where the vis- 
cous and the buffer layers cannot be resoled and should 
be parameterized. To attain these ends we need to ob- 
tain analytical solutions (numerical solution would be 
of no use), which calls for simplification of the govern- 
ing equations. Accordingly, our strategy is a pragmatic, 
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task-dependent simplification and restrictions. In par- 
ticular we concentrate on descriptions of the profile of 
mean flow and the statistics of turbulence on the level 
of simultaneous, one-point, second-order velocity corre- 
lation functions. In other words, the objects that we are 
after are the entries of the Reynolds-stress tensor as a 
function of the distance from the wall. The model will 
be presented for plain geometry; this geometry is rele- 
vant for a wide variety of turbulent flows, like channel 
and plain Couette flows, fluid flows over inclined planes 
under gravity (modelling river flows), atmospheric tur- 
bulent boundary layers over flat planes and, in the limit 
of large Reynolds numbers, many other turbulent flows, 
including pipe, circular Couette flows, etc. 

Suggested in this paper phcnomenological theory of 
wall-bounded flows is based on standard ideas [5] ; never- 
theless we develop the theory slightly further than any- 
thing that exists currently in the literature. In our study 
we will stress analytical tractability; in other words, 
we will introduce approximations in order to achieve a 
model whose properties and predictions can be under- 
stood without resort to numerical calculations. Never- 
theless we will show that the model appears very de- 
pendable in the sense that its predictions check very well 
in comparison to direct numerical simulations, including 
some rather non-trivial predictions that are corroborated 
only by very recent simulations and experiments (which 
only now reach the sufficient accuracy and high Reynolds 
numbers). 

We should notice, that considering the mean velocity 
and the second-order statistics in these (neutrally strat- 
ified) flows, we neglect some mechanisms and features 
although present but not essential in the problem under 
consideration. However, proceeding further, in particu- 
lar, to account for the density or temperature stratifica- 
tion (out of the scope of this paper), we quite probably 
will be forced to rule out of some simplifications accept- 
able in the first task and to account, for example, for a 
spacial energy flux and even for coherent structures. 

In Sec. I we formulate a model which is a version of the 
"algebraic Reynolds-stress model" [5]. In Sec. I A we in- 
troduce notations and recall the equations describing the 
mechanical balance; in Sec. IB we state the assumptions 
and detail the approximations used in the context of the 
balance equations for the components of the Reynolds 
stress tensor Wij . The result of these considerations is a 
set of 5 equations for the mean shear 5* and Wij which is 
described in Sec. IC. For actual calculations this set of 
equations is still too rich since it contains 12 adjustable 
parameter. Eight of these parameters control the non- 
linear behavior of the system in the outer layer and four 
additional parameters govern the energy dissipation in 
the viscous sub- layer. Clearly, further reduction of the 
model is called for. This is accomplished in Sec. II. 
First, in Sec. II A we consider the full 12 parametrical 
solution of the model, and present a comparison with ex- 
perimental observations in Sec. II A 4. This comparison 
indicates that an adequate description of the entire tur- 



bulent boundary layer phenomenology can be achieved 
with only four parameters instead of twelve. We refer 
to the four-parameter model as the "minimal model" . 
In Sec. II B we reap the benefit of the minimal model: 
we find simple and physically transparent Eqs. (2.15) for 
the profiles of the Reynolds stress tensor Wij (y) and the 
mean shear S(y) (y is the distance from the wall) in terms 
of the root-mean-square turbulent velocity v = y/Wu. 
Unfortunately, the equation for the v(y) profile is quite 
cumbersome and cannot be solved analytically. Never- 
theless we employ an effective iteration procedure that 
allows reaching highly accurate solutions with one or at 
most two iteration steps. 

Section III is devoted to a comparison of the pre- 
dictions of the minimal model with results of experi- 
ments and direct numerical simulations . In particular, in 
Sec. Ill B we show that the model describes the mean ve- 
locity profile in a channel flow with 1%-accuracy almost 
everywhere. Only in the core the model fails to describe 
so-called "velocity defect" (the upward deviation from 
the log-law) which is observed near the mid-channel (of 
about 5-6 units in V + , independent on Reynolds num- 
ber). For our purposes this mismatch in not essential. 
In Sec. Ill C we show that the minimal model provides a 
good qualitative description of kinetic energy profile, in- 
cluding position, amplitude and width of the peak of the 
kinetic energy in the buffer layer. In Sec. Ill D we show 
that with the same set of four parameters the model of- 
fers also a good qualitative description of the Reynolds 
stress profiles and the profiles of "partial" kinetic energies 
(in the streamwise, wall-normal and cross-stream direc- 
tions) almost in the entire channel. The final Sec. IV 
presents a short summary of our results, including a dis- 
cussion of the limitations of the minimal model. Possible 
improvements of the suggested model will have to start 
by addressing these limitations. 

I. FORMULATION OF THE MODEL 

Our starting point is the standard Reynolds decom- 
posion [5] of the fluid velocity U(r,t) into its average 
(over time) V and the fluctuating components u. In 
wall-bounded planar geometry the mean velocity is ori- 
ented in the (stream-wise) x direction, depending on the 
vertical (wall-normal) coordinate y only: 

U(r, t) = V(y) + u{r, t) , V(y) = (U(r, t)) = x V(y) . 

(1-1) 

The mean velocity and the fluctuating parts are used to 
construct the objects of the theory which are the compo- 
nents of the Reynolds stress tensor W(y) and the mean 
shear: 

Wij^^^Uj) , S(y)=^-. (1.2) 

We note that in previous applications [6-14] we have em- 
ployed a model in which only the trace of W(y) and its xy 
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component were kept in a simplified description. For the 
present purposes we consider all the component of this 
tensor, paying a price of having more equations to bal- 
ance, but reaping the benefit of a significantly improved 
phenomenology. We discuss now the equations relating 
these variables to each other. 



A. Equation for the mechanical balance 

The first equation relates the Reynolds stress W xy to 
the mean shear; it describes the balance of the flux of 
mechanical momentum, it follows as an exact result from 
Navier-Stokes equations and has the familiar form: 

-W xv (y) + v S{y) = P(y) . (1.3a) 

W xy on the left hand side (LHS) is the turbulent (re- 
versible) contribution to the momentum flux whereas 
vqS(d) is the viscous (dissipative) contribution to the 
momentum flux. The RHS is the momentum flux, which 
may have different origin. For example, in a channel or 
pipe flow P{y) is generated by the pressure head. In the 
channel flow with the pressure gradient p' = —dp/dx, 

P(y)=p'(L-y), (1.3b) 

where L is the half width of the channel. In the pipe 
flow P(y) is given by the same Eq. (1.3b) with L being 
a half of the pipe radius. In a water flow over incline in 
the gravity field p' should be replaced by gsina, where 
g is the gravity acceleration and a is the inclination an- 
gle. For Re a S> 1, near the wall one can neglect the y 
dependence of P(y), replacing P(y) by its value at the 
wall: P(y) => P = P(0). Here the so-called "wall-based" 
Reynolds number Re a is defined by: 

RCJ . him . (1 .4) 



B. Balance of the Reynolds tensor 

The next set of equations relates the various compo- 
nents of the Reynolds tensor, Wij(y) defined by Eq. (1.2). 
In contrast to Eq. (1.3a) this set of equations is only par- 
tially exact. We need to model some of the terms, as 
explained below. We start from the Navier-Stokes equa- 
tions and write the following set of equations: 

^jf- + + J ij = S(W iy 6 jx + W jy S tx ) . (1.5a) 

The RHS of these equations is exact, describing the pro- 
duction term in the equations for Wi X = W X i which is 
caused by the existence of a mean shear. On the LHS of 
Eq. (1.5a) 

/ dui duj \ . „, . 

e « = 2 "°Wsr) (L5b > 



is the exact term presenting the viscous energy dissipa- 
tion. The problem is that involve new object, which 
requires evaluation via Wij. This can be easily done in 
regions where the velocity field is rather smooth, and in 
particular in the viscous sub-layer, the velocity gradient 
exists and thus the spatial derivatives in Eq. (1.5b) are 
estimated using a characteristic length which is the dis- 
tance from the wall y. In order to write equalities we 
employ the dimensionless constants ~ 1: 

ei^e^=7^, 7™^(^) 2 - (l-5c) 

In general the constants ay are different for every 

In the buffer sublayer and in the log-layer the energy 
cascades down the scales until it dissipates at the Kol- 
mogorov (inner) scale that is much smaller than the dis- 
tance y from the wall. Therefore the main contribution 
to the dissipation from all scales smaller than y is due 
to the energy flux, i.e. has a nonlinear character. Due 
to the asymptotical isotropy of fine-scale turbulence, the 
nonlinear contribution should be diagonal in i, j (see, 
e.g. [5]): 

W 

e ij ^ef j = 1 — S ij , W^Tt{W}, (1.6a) 

where prefactor i is introduced to simplify equations be- 
low. The characteristic "nonlinear flux frequency" j nl , 
can be estimated using a standard Kolmogorov-41 di- 
mensional analysis: 

i(y) = -Vw(yj, , (i-6b) 
y 

again with some constants b <~ 1. The "outer scale" of 
turbulence is estimated in Eq. (1.6b) by the only available 
characteristic length, y, the distance to the wall. 

As one sees from Eq. (1.6), the dissipation of partic- 
ular component of the Reynolds-stress tensor, say W xx , 
depends not only on W xx itself, but also on other com- 
ponents, W yy and W zz in our case. It means that ey, 
given by Eq. (1.6a) leads, in the framework of Eq. (1.5) 
not only to the dissipation of total energy, but also to 
its redistribution between different components of Wu. 
In order to separate these effects we divide ey into two 
parts as follows: 

4 1 = . (l-7b) 

4 2 = -7(^-%y)% • (1.7c) 

Clearly, e^' 1 describes the damping of each component 
Wu separately, without changing of their ratios, while the 
traceless part, e^-' 1 , does not contribute to the dissipation 
of total energy and leads only to redistribution of energy 
between components of the Reynolds-stress tensor. This 
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contribution we will include into the "return to isotropy" 
term (1.14b), that will be discussed below. 
Actually, we presented as the sum 



^dis 



al,2 



(1.8) 



in which for the total energy dissipation is responsible 
only first term in the RHS. In the buffer layer both con- 
tributions to e^ ls , the viscous dissipation, 
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Hj ' 



nonlinear one, e?-' are important. Their relative role de- 
pending on the turbulent statistics. We will employ two 
simple interpolation formulas which lead to two versions 
of the minimal model: 



r u (y) 
r ij (y) 



i(y) k 



'YiiM S (y) + 7 2 (y)k 



(1.9) 

"sum" , (1.10) 
"root" .(1.11) 



The versions of the resulting model will be referred to as 
the "sum" and the "root" versions correspondingly. A- 
priori there is no reason to prefer one or the other, and 
we leave the choice for later, after the comparisons with 
the data. 

The term in Eq. (1.5a) is caused by the pressure- 
strain correlations: 



hj — 



1 

Po 



dxj 



dxi 



(1.12) 



and is known in the literature as the "Return to Isotropy" 
[5]. Due to incompressibility constraint Jjj is traceless 
tensor and therefore does not contribute to the total en- 
ergy balance, leading only to redistribution of partial 
kinetic energy between different vectorial components. 
Also, this term does not exist in isotropic turbulence 



where W;. 
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We adopt the simplest linear Rota 



approximation for the "Return to Isotropy" term [5], 
using yet another different characteristic frequencies 7 ^ , 
estimated as follows: 



liAy) 



7y(3Wtf-5«W) 



(1.13a) 
1 . (1.13b) 



Thus, representation (1.14b) involves only two free pa- 
rameters bd and b. 

Equations (1.5) for with ij = xx, yy, zz and 

xy involve 7 constants ay, 6, &d and b. Our goal is to 
formulate the simplest possible model, with a minimal 
number of adjustable constants. The strategy will be now 
to use experimental and simulational data, coupled with 
reasonable physical considerations, to reduce the number 
of parameters to four, each of which being responsible for 
a separate fragment of the underlying physics. 

We should stress that we neglect in Eq. (1.5a) the 
spatial energy transport term e tr , caused by the tripplc- 
velocity correlations, pressure- velocity correlations and 
by the viscosity [5]. In the high Rca limit the density 
of turbulent kinetic energy becomes space independent 
in the log-law region. Accordingly, the spatial transport 
term is very small in that log-law region. More detailed 
analysis, see, e.g. Fig. 3 in Ref. [15], shows that even for a 
relatively small ReA in the log-law turbulent region this 
term is small with respect to the energy transfer term 
from scale to scale which is represented by jWij in the 
equations above. On the other hand, in the viscous sub- 
layer the mean velocity is determined by the viscous term 
and thus the influence of the spatial energy transfer term 
can be again neglected. To keep the model simple we will 
neglect e tr term also in the buffer layer where it is of the 
same order as the other terms of the model. The reason 
for this simplification, which evidently will cause some 
trouble in the buffer layer, is that the energy balance 
equations used below become local in space. This is a 
great advantage of the model, allowing us to advance an- 
alytically to obtain a very transparent phenomenology of 
wall-bounded turbulence. It was already demonstrated 
in Ref. [12] that the simple description (1.10) gives a 
uniformly reasonable description of the rate of the en- 
ergy dissipation in the entire boundary layer. Here we 
improve this description further, effectively accounting 
for the energy transfer term in the balance equation by 
an appropriate decrease in the viscous layer parameters 



One sees that hj has precisely the same structure as e"j ' 2 , 
introduced by Eq. (1.7c). Therefore it is convenient to 
treat these contributions together, introducing 



7 J -A- n1,2 

J-ij — J-ij T Hj , 

4- = % (3Wij - SijW) , 



(1.14a) 
(1.14b) 



y 



hj = hj - - (% . (1.14c) 



Recall that tensor 1^ must have zero trace for any val- 
ues of Wij. This is possible only if b xx = b yy = b zz = 
fed , b xx = 6, and consequently 



lyy = Izz = Id , Ixy = 7 • 



(1.15) 



C. Summary of the two versions of the model 

For the sake of further analysis we present the model 
with the final notation: 



-W xy {y) + v S{y) 


= P(y), 


(1.16a) 


\Pxx + &rd]w xx 


= ^ d W - 2SW xy , 


(1.16b) 


[Tyy +3j d ]Wyy 


= IdW, 


(1.16c) 




= IdW, 


(1.16d) 


[T xy + 3j]W xy 


= SWyy, 


(1.16c) 



In the traditional theory of wall-bounded turbulence one 
employs the "wall units" u Ti t and l T for the velocity, 
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time and length [5] which for a fluid of density p are: 



(1.17) 



Using these scales one defines the wall-normalized dimen- 
sionless objects 

VHy) = ^, vl { y)^-M, etc., 



V 

s+ 



V 



St. 



etc. 



(1.18) 



In our model we can use the property of locality in space 
to introduce "local units" : 



u T {y) 



, r{y) 



P{y) 



, tr(y) 



VppW)' 

(1-19) 

similar to traditional wall units Eq. (1.17), but with the 
replacement Pq — > P(y), and "locally normalized" di- 
mensionless objects, analogous to Eq. (1.18): 



£ T (y) ' " My) 

w J3 {y) 

u%(y) 



(1.20) 



S*(y*) ee Sr(y), W^tf) 



etc. . 



Then the dimensionless version of Eq. (1.16a) reduces to 



-wL + s* = 1 



xy 



(r^ + 3 7dX 
(rt w + 37*)wt B 

Introducing = VW$ we can write: 



7^-2^4, 



(1.21a) 
(1.21b) 
(1.21c) 
(1.21d) 
(1.21e) 



it 



rt. 



7 d 



(^r) 2 + (5i 7 , for the sum model , (1.22a) 

y x y % 



1 a 4 . b 2 vt 

J^y + ' f ° r th ° r00t modcl (!-22b) 



bdv t 



— , for both versions .(1.22c) 



II. ANALYSIS OF THE MODEL 

A. Solution of the 7-parameters version of the 
model 

1. Solutions in the viscous sub-layer 



W zz and W xy (with coefficients that are functions of 
W). They can have a trivial solution W = for which 
Eq. (1.16a) gives 

S = Po o j V = yPo I v a Wij = , => Laminal layer. 

(2.1) 

The complete absence of turbulent activity in the viscous 
layer in our modcl is a consequence of leaving out the 
energy transport in physical space. 



2. Analysis of turbulent solution 

Equations (1.21b - 1.21e) have non-trivial "turbulent 
solution" with W t ^ 0: 



W x = 



r* 



7 d 



"yy 



xy 



: jL_ w t _^*M 

" z. 



(2.2a) 



(rI, + 37*)(rL + 373 



♦ , (2.2b) 

rL + 3^ 

(2.2c) 



if its determinant A vanishes. The solvability condition 
A = gives: 



(Si) 2 = 



xy 1 / xy 



yy y- zz + 3 "f z 



rLrLrt (2.2d) 



xx yy zz 



27 

xo^t/'rl pi _i_ pi pi _i_ pi 

'd V a:a: yy 1 a; a; 2:2: 1 yy x 2:2: ,/ 

+3 (7^) 2 (rL + r^ + rJ 



Substitution M^J and S$ in Eq. (1.21a) gives a closed 

equation for the function W$(y$), (or for = V\V$). 
To present the resulting Eq. in explicit form, introduce 

A(v*) = S^)^-^,/^^, (2.3) 

Using Eqs. (2. 2d) and (2.2c) we find: 



AV) 



3r 



2r d (R zz + 3r d 

~f~2 7*d (RxxRyy ~l~ RxxRzz 
Pyy Pzz) 



RxxRyyRzz (2-4) 



RyyRzzj 



+3r 2 jR x 



TO + O Td J 

Now Eq. (1.21a) can be presented as 



(2.5) 



The four Eqs. (1.16b) - (1.16e) can be considered as 
a homogeneous "linear" set of equations for W xx , W yy , 



Together with Eqs. (2.2) this provides the full solution 
of Eqs. (1.21). 
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3. Outer layer, j/* > 50 

In the outer layer, far away from the wall, all the vis- 
cous terms in Eqs. (1.21) can be neglected. In this case 
Eqs. (1-22) for both the sum and the root minimal models 
give: 



-4 



7^, 



(2.6) 



and Eqs. (2.2b-2.2b) and the solution (2.15) simplifies 
drastically: 



W t = W X = 

VV zz 



wHd_ 

b + 3b d 



(2.7) 



By analyzing results of experiments and numerical sim- 
ulations (as discussed in Sec. Ill) we found that in the 
outer layer W x x = \W X , and W* y = W ZZ = \W % . The 
model reproduces these findings if we choose: 



b = b d . 



(2.8) 



Using this relation and solution of (1.21a): W x y = —1, 
in the rest of Eqs. (1.21), one finds 



W x = 




1 

K 



1/4 



(3 6) 3/4 . (2.9) 



Here k is nothing but the von-Karman constant, that 
determines the slope of the logarithmic mean velocity 
profile in the log-law turbulent region: 

V+(y+) = K- 1 lny+ + C, forz+>30, (2.10) 
k w 0.436, C w 6.13 . 

The experimental value of k and the intercept C, were 
taken from [18]. Using the simulations result W + s=s 6.85 
of Rcf. [15] which is reproduced in Fig. 6, we find 



0.256. 



0.500 . 



(2.11) 



4- Reduction of the number parameters: the minimal model 

The parameters ay are responsible for the difference 
between the energy dissipation and the energy transfer in 
the viscous sub-layer. To further simplify the model we 
reduce the number of independent parameters ay from 



4 (a 



xx i ^yyj ^zz 



Among various possibilities (including a xx — a zz = a, 

®"yy — ^? Qxy 

= (a + a)/2) we choose a parametrization 
similar to the situation with the outer layer parameters: 

an — a , a xy = a . (2-12) 

The analytical solution given in next Sec. II B simplifies 
considerably with the 4-parameter version of the model. 
A further simplification a = a could be considered, but 
we rule it out since it yields a monotonic dependence 



of the turbulent kinetic energy W + {y + ) with y + , while 
experimentally there is a pronounced peak of W + {y + ) in 
the buffer sub-layer, see Sec. IIIC. We thus consider the 
four-parameter model as the "minimal model" (MM). 
Below we will use mostly the following set of constants: 



a = 1.0, 5=10.67, sum-MM, (2.13a) 

a= 1.0, 5=12.95, root-MM, (2.13b) 

b = 0.256, 6 = 0.500, both MMs. (2.13c) 

This choice is based on the analysis of the simulational 
and experimental data presented in Sec. III. 

Notice that eliminating b from Eqs. (2.9) (valid for 
both sum- and root models) one gets: 



v x = 12nb w 6 k . 



(2.14) 



With the simulational values k — 0.436 and W X oo = 6.85 
this relationship is valid with a precision that is better 
than 1%. 



B. Analysis of the Minimal Models 

1. a,a-parametrization of the general solution 

With the minimal parametrization, given by 

Eq. (2.12), the solution (2.2) takes on a simpler 
form: 

W* = W\ z = ^-W i , W* x = W x , (2.15a) 

yy 4u4 > xx 2v 4 y ' 



W* = - 



xy 



W x I bv^vi ^ _ 



2 V 66 



V3V4 



y 



LJ^p-i .(2.15b) 



Here we introduced the following short-hand notations 



Vj for the sum-MM: 



Vi = v x 



byt> 



v 2 = + 



2byt 



(2.16a) 



^3 = v l H — - — , V4 = v l + 
3byt 



4byi ' 



For the root-MM instead of Eq. (2.16a) we take: 



and a xy ) to two, denoted as a and a. Vl 



2 + ^— 

(by*)' 



lvt 2 + 



, v 2 



v 3 = ^ + . {b-bfvt 2 + " 



vi +v x vi + 3v x 

"~ 4 ' 

(2.16b) 



(Uytf 



With the minimal parametrization Eq. (2.5) takes a 
very simple explicit form: 



t 2 , I2bv % d j24br 3 r 4 



w + + 



r 3 r 4 



bri 



rj = vj/vi . (2.17a) 
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This form of equation for uj. serves below as a starting 
point for an approximate (iterative) analytical solution. 
One can also seek an exact solution by numerical meth- 
ods; to this aim it is better to use the following form of 
the same equation: 



A much better approximation for v*(y*) (denoted as 
v\) is obtained using in Eq. (2.20) a -independent r^ = 
r j( v io) instead of rj = 1: 



b_ 

24 6 



F(v*, y*) = -^t)'V + A 3 »4 



6 t 2 
— v* v\ 













4 = \ 


/ 24 6 r 3 . r 4i0 , 


( 6 6r 3 . r 4 . \ 


1 


V br lfi + 


[ y* ) 



6bb 22 
H — v 1 v 3 v i = . 



(2.17b) 



Equation (2.17b) has seven roots for the sum-MM, 
(and 27 for the root-MM) but only two of them, denoted 
as v± are real and positive for large enough y* . These two 
roots approach each other upon decreasing the distance 
from the wall. At some value of y* these roots merge: 



6 6r 3 . r 4 . 

v* ' 

(2.22a) 

Clearly, this iterative procedure can be prolonged further 
and one can find the velocity v$ at the n + 1 iteration 
step, vl l+1 (y$), using the relations r^ n = Tj{v\), found 
with the velocity v\ of the previous, n-th step: 



4+i 







/ 24 b r 3;n r 4i „ 


( 6 6r 3i „r 4 .„\ 


V br l,n 


I y* ) 



6 6r 3 „r 4i „ 



y 



t 



v+iyvs) = vi(y vs ) = w» < 



(2.18) 



The values y va and v* as functions of the problem param- 
eters follow from the polynomial (2.17b): 



F(v t ,y t ) = o, 



dF(vt,yt 
dvt 



= 



(2.19) 



For y* < y vs there are no physical (positive definite) 
solutions of Eq. (2.17b). This is a laminar region that 
was discussed before as the viscous sub-layer. In Tab. II 
we present the corresponding values of y vs and for 
b = 0.256 , b = 0.5 and various pairs of a , a. 



2. Iterative solution of Eq. (2.17a) for rms turbulent 
velocity v*(y*) 



To develop further analytic insight we employ an itera- 
tive procedure to find an approximate solution for u+(y) 
for all y^ > y vs . For this goal we forget for a moment that 



(2.22b) 

The numerical verification of the iteration procedure is 
given in Appendix A. The conclusion is that already 
the first few iterations are sufficiently accurate for all 
practical purposes: often one can use the first iteration 
and occasionally the second one. 

Remarkably, the first iteration can be formulated di- 
rectly in terms of the basic Eqs. (1.21) by replacing 
the turbulent velocity profile v*(y*) in Eqs. (1.22) for 
and ■yjj by its asymptotic value in log-law region 

1/4 



>t = 



(24 6/6 



3. Iterative solution for the mean velocity and Reynolds 
tensor profiles 

Consider first the resulting plots for the mean veloc- 
ity profile, V*(y^), computed with the help of turbulent 
velocity v\ at the n-th iteration step: 



■ all y* > y vs . for tms goal we lorget tor a moment tnat rv 

depends on v x , and consider Eq. (2.17a) as a quadratic = 2/vs + / S*(0<Z£, y X > y vs ■ (2.23) 

...j.: ja - 



equation with a positive solution: 



' 24 6 r 3 r 4 j 6 6r 3 r 4 



6 6 r 3 r 4 



bri 



y 



(2.20) 



However rj does depend on v*. For example, for the 
sum-MM: 



ri(^) = 1 + 



1 + 



byivt ' 
3 6y : M 



r 2 (v t ) = 1 + 
r 4 (w*) 



1 



2 6yM 
a 2 
4 6yM 



, (2.21) 



Nevertheless, for very large y* all rj — > 1 and an asymp- 
totic solution of (2.20) reproduces the asymptotic value 
of v x = ?4 = (246/6) 1 /4 j gi vcn by Eq. (2.9). 



Here S$(x) denotes S$(x), given by Eqs. (2.15), with = 
u*. Figure 1 displays plots of V^(y^) for n = 1, 2, 3, 4 
and the "exact" (numerical) result V*(y*). All the plots 
almost coincide within the linewidth. This means that 
for the purpose of computing V$(y$) one can use the 
first approximation v\ given by Eq. (2.22a) instead of 
the exact solution v*. 

Next we present in Fig. 2 log-plots for the trace of the 
Reynolds-stress tensor W\ (y* ) , (computed with the n-th 
iteration step for n = 1, 2, 3, 4) together with the "exact" 
numerical solution W*(y*) for the sum-MM. Evidently, 
the iterative procedure for the kinetic energy does not 
converge as rapidly as for the mean velocity profile: one 
can distinguish the plots of W} (y*), W\ (y*) and W 3 (y*); 
the plots of W 3 (y*) and wj(y^) coincide within the line 
width. Nevertheless, for y* > 5 (i.e. in the buffer layer 
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by C. Casciola [16], and to recent laboratory experiments 
in a vertical water tunnel by A. Angrawal, L. Djenidi and 
R.A. Antonia [17]. The choice of the outer layer parame- 
ters bij and bij is based on our analysis of the anisotropy 
in the log-low region, presented in Sec. Ill A. The relation 
between the viscous layer parameters a vs. a is based on 
the comparison between the DNS and the model mean 
velocity profiles, presented in Sec. IIIB. The final choice 
of a and a is motivated by the DNS data for the kinetic 
profile which is compared with the model prediction in 
Sec. IIIC. Section HID is devoted to the comparison of 
the model results with the DNS profiles of the Reynolds 
stress Wl y and partial kinetic energies W* x , W* y , W\ z . 



FIG. 1: Log-plots of the exact solution for mean velocity pro- 
file V*(y*) and approximate profiles V^(j/*), computed on the 
n-th iteration step for n = 1,2,3,4 for sum-MM with con- 



++ c 




FIG. 2: Log-plots of the iterative W„, n = 1 — 4 and "exact" 
numerical solutions (thick solid line) for trace of the Reynolds 
tensor in sum-MM with constants, taken from Eq. (2.11). 
Plots for W3 and W\ coincide within the line width. 



and in the outer region) alreadly the first iterative so- 
lution provides a very reasonable approximation to the 
exact solution for the kinetic energy profile. 



III. ANALYSIS OF THE NUMERICAL AND 
EXPERIMENTAL DATA AND COMPARISON 
WITH THE MODEL PREDICTION 

In this Section we analyze and compare the predic- 
tions of the minimal models to results of experiments and 
comprehensive direct numerical simulations of high Re 
channel flows. We refer to results that were made avail- 
able in the public domain by R. G. Moser, J. Kim, and 
N. N. Mansour [15], to Large-Eddy-Simulation performed 



A. Anisotropy of the log-layer: Relative partial 
kinetic energies R xx , R yy , Rzz in the outer layer 

The anisotropy of turbulent boundary layer can be 
characterized by the dimcnsionless ratios 



Ru(y + ) 



W(y+ 



W+(y+) 



(3.1) 



This anisotropy plays an important role in various phe- 
nomena and was a subject of experimental and theoret- 
ical concern for many decades, see, e.g. [1, 5]. Never- 
theless, up to now the dispersion of results on the sub- 
ject appears quite large. There is a widely spread opin- 
ion, based on atmospheric measurements, that the wall- 
normal turbulent fluctuations W yy are much smaller than 
the other ones. For example, Monin and Yaglom [1] re- 
ported that for a neutrally stratified log-boundary layer 
R xx = 54%, R yy = 6% and R zz = 40%. This contra- 
dicts recent DNS results for ReA=590 which are avail- 
able in Ref. [15], as shown in Fig. 3. Note that there 
is a region about 100 < y + < |Re> where the plots of 
Ru{y + ) are nearly horizontal, as expected in the log-law 
region. From these plots we can conclude that is this re- 
gion R xx ss 53% which is close to the value 54%, stated 
in [1]. Nevertheless, the DNS data for R yy and 1Z ZZ are 
completely different. From Fig. 3 one gets TZ yy ~ 22% 
and 1Z ZZ « 27%. Thus lZ yy can be considered roughly 
equal to 1Z ZZ . We should mention here that various mod- 
els of turbulent boundary layers give lZ yy — 1Z ZZ in the 
asymptotic log-law region. We propose that the differ- 
ence between lZ yy and 1Z ZZ which is observed in Fig. 3 
is due to the effect of the energy transfer. This effect 
practically vanishes in the asymptotic limit Re a — > 00, 
but is still present at values of ReA which are available 
in DNS [15]. Indeed, for both values of ReA shown in 
Fig. 3, Wyy = W zz in the center of the channel, where 
the energy flux vanishes by symmetry. Clearly, there is 
no energy flux also in space homogeneous cases, for ex- 
ample for a constant shear flow, in which, according to 
the model, one should expect W yy — W zz in the entire 
space. 

Our expectation that W yy = W zz , which is based on 
symmetry considerations, is confirmed by the Large Eddy 



Simulation (LES) of the constant shear flow [16]. As one 
sees in Fig. 4 in this flow 1Z XX w 0.46, while W yy w W zz w 
0.27. 

As stated, for sufficiently large values of Re>, the en- 
ergy transfer terms should almost vanish in the log-law 
region and, according to our model, one can expect in 
that region W yy = W zz also in the channel flow. This 
viewpoint was confirmed in the aforementioned labora- 
tory experiment [17] in a vertical water channel with 
Rca = 1000, reproduced in Fig. 5. The experimental val- 
ues of Ttxxi T^-yy an d TZ-zz m the log- law turbulent region 
are in the excellent quantitative agreement with the val- 
ues 1Z XX = 0.5 and lZ yy = 1Z ZZ = 0.25 shown in Figs. 3 
and 5 by horizontal dashed lines. 

Table I summaries the DNS, LES and experimental 
values of the relative kinetic energies in comparison with 
the model expectation. The conclusion is that, in contra- 
diction with the old and still wide spread viewpoint [1] 
that the wall-normal turbulent activity is strongly sup- 
pressed, lZyy < 0.1, the turbulent kinetic energy in the 
log-law region is distributed in a very simple manner: the 
stream-wise component contains a half of total energy, 
T^-xx = h and the rest is equally distributed between the 
wall-normal and span- wise components: 7Z yy — 7Z ZZ = \- 
As shown in Sec. II A 4, this very simple energy distribu- 
tion is predicted by the minimal model if one assumes 
that the characteristic nonlinear times scales in the en- 
ergy transfer term and in the return-to- isotropy term are 
identical. 

Thus anisotropy predicted by our minimal model agree 
reasonably accurately with those obtained from DNS, 
LES and vertical water channel. However, we must ad- 
mit that all these are not yet the nature. Indeed DNS, 
LES and lab experiments, done at relatively modest ReA 
impose limits on the low-frequency intervals in the spec- 
tra of the streamwise and the transverse velocity compo- 
nents, because of side walls or periodic boundary condi- 
tions. In other words, it must not be ruled out that DNS, 
as well as lab experiments, cat off the largest-scale ejec- 
tions, observed in the atmosphere in the form of coher- 
ent structures, which pump additional streamwise- and 
transverse-velocity energies into the log-layer. Be it as 
it may, we leave a detailed discussion of the above prob- 
lem and geophysical applications of our theory for fur- 
ther work. At the present stage, following our strategy 
of "pragmatic, task-dependent simplification", we con- 
sider our minimal model as definitely relevant to flows in 
channels. Its extension to geophysical (atmospheric and 
oceanic) boundary layers needs further efforts. 



B. Mean velocity profile in channel flows 

To compute the mean velocity profile V + (y + ) in our 
approach we need to connect first S + (y + ) with S*(y*). 




+ 

y 



FIG. 3: DNS profiles of the relative kinetic energies in the 
stream- wise, wall-normal and span- wise directions, R xx , R yy 
and R zz respectively. Solid lines: ReA = 590, dot-dashed 
lines: ReA = 395. Horizontal dashed lines show levels 0.5 and 
0.25. 
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FIG. 4: Relative components of the Reynolds stress tensor 
Wu, Eq. (3.1), for for constant-shear flow. Results of the 
LES [16]. Lines serve only to guide the eye. 



TABLE I: Asymptotic values of the relative kinetic ener- 
gies Ha in the log-law region y* > 200 taken from DNS 
(ReA = 395, 590) [15], LES [16] and experiment in a water 
channel with ReA = 1000 [17]. The last column presents the 
predictions of the minimal model. 





DNS [15] 


LES [16] 


Water channel [17] 


Model 


XX 


« 0.53 


« 0.46 


0.50 ±0.01 


0.50 


yy 


« 0.22 


« 0.27 


0.25 ±0.02 


0.25 


zz 


« 0.27 


« 0.27 


0.25 ±0.02 


0.25 
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FIG. 5: Experimental profiles of the relative kinetic energies 
lZa in a vertical water channel with ReA = 1000 according 
to Ref. [17]. The solid lines serve to guide the eye. The 
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V + {y + ) for ReA = 590 and the parameters (2.11) are 
shown in Fig. 6 as a dashed line for the sum-MM and as a 
dot-dashed line for the root-MM. The DNS profile of [15] 
for the same ReA is shown as a solid line. There is no 
significant difference (less than 1%) between these plots 
in the viscous sublayer, buffer and outer layers, where 
y + % 300 i.e. in about 50% of the channel half-width 
L + = Re a = 590. This robustness of the mean velocity 
profile V + (y + ) is a consequence of the fact that V + (y + ) 
is an integral of the mean shear S + which is described 
very well both in the viscous and the outer layers. 

Notice that our model does not describe the upward 
deviation from the log-low which is observed near the 
mid-channel (of about 5-6 units in V + , independent on 
Reynolds number). We consider this minor disagreement 
as an acceptable price for the simplicity of the minimal 
model which neglects the energy transport term toward 
the centerline of the channel. This transport is the only 
reason for some turbulent activity near the centerline 
where both the Reynolds stress W xy and the mean shear 
S vanish due to symmetry. Just at the center line the 
source term in our energy equation, —2SW xy , is zero, 
and the missing energy transport term is felt. 

The plots in Fig. 6 have a reasonably straight logarith- 
mic region from y + ss 20 to y + s» 200. On the other 
hand, the Reynolds stress profile at the same ReA = 590 
shown in Fig. 8, has no fiat region at all. Such a 
flat region is expected in the true asymptotic regime of 
ReA — > oo, where W + = — 1. Therefore if one plots the 
model profiles V + at different ReA and fits them by log- 
linear profiles (2.10) one can get a ReA-dependence of 
the "effective" intercept in the von-Karman log-law. We 
think that this explains, why measured value of the log- 
low intercept can depend on the Reynolds number and 
on the flow geometry (channel vs. pipe): both in DNS 
and in physical experiments one usually does not reach 
high enough values of ReA ■ 



FIG. 6: Mean velocity profiles V + (y + ): The solid blue line - 
DNS data [15] for Re T = 590, the dashed green line - analyti- 
cal profile V + (y + ) for the sum-MM, and dash-dotted red line 
for the root-MM using the parameters Eq. (2.11). The rela- 
tive deviation of the analytical predictions from the DNS data 
are about 1%, smaller for the root-MM than for the sum-MM. 



According to the definitions (1.18) - (1.20): 



S + (y + ) = 




where 

ReA = L/£ T , (3.3) 

and y + > y vs . For y + < y vs we can take S + (y + ) = 1 and 
integrate the resulting shear over the distance to the wall 
with no-slip boundary condition. The resulting profiles 



C. Profiles of the total kinetic energy density and 
the choice of the pair a, a 

The quality of the profiles V + (y + ) calls for a bit more 
thinking. In fact, one find that the minimal model pro- 
duces practically the same profiles V + (y + ) not only for 
the parameters (2.13) but for a wide choice of the pairs 
a, a, for example for a — 2 and a — 8.6. Actually, for 
any < a < 4 one can find a value of a that gives a 
mean velocity profile in good agreement with Fig. 6. In 
other words, in the (a, a)-plane there exist a long narrow 
corridor that produces a good quantitative description 
of V + (y + ). Within this corridor there exists a line that 
provides a "best fit" of V + (y + ), minimizing the mean 
square deviation 6V + 

5V+ = ^/([V+ M (y + )-V^ s (y+W) (3.4) 

of the model prediction V+ M (y + ) from the DNS profile 
Vj NS (y + ) in the inner region y + < 140. Some of the best 



TABLE II: Optimal values of a for a given value of a, that 
JV+ofEq. (3.4). For optimal pairs (a, a) in the sum and the 
root versions of the minimal model (denoted as an( i 
correspondingly) we present the values y va of t/*, separating 
the viscous solution with v* — from turbulent regime with 
v* = v*(y VB + 0). The last two columns present the j/* po- 
sition, j/max, of the maximum of the kinetic energy and the 
corresponding values of W,* ax = W*(y ma x). 



MM 


a 


a 


8V t 




«* 


J/max 




E 


0.1 


10.4 


0.29 


1.5 


0.017 


21 


8.33 


v~ 


0.1 


12.0 


0.25 


l.i 


0.008 


18 


8.32 


E 


0.25 


10.9 


0.25 


2.4 


0.061 


22 


8.48 


V~ 


0.25 


12.1 


0.25 


2.1 


0.016 


19 


8.36 


E 


0.5 


11.1 


0.25 


3.4 


0.159 


22 


8.47 




0.5 


12.6 


0.17 


3.4 


0.247 


19 


8.50 


E 


1.0 


10.7 


0.22 


4.8 


0.401 


24.6 


8.24 


V~ 


1.0 


12.9 


0.15 


4.9 


0.234 


19.2 


8.59 


E 


1.5 


9.7 


0.20 


5.7 


0.600 


27 


7.76 


V~ 


1.5 


12.9 


0.16 


4.8 


0.743 


19 


8.58 


E 


2.0 


8.6 


0.19 


6.4 


0.783 


27 


7.26 




2.0 


11.8 


0.17 


6.7 


0.743 


20 


8.13 


E 


4.0 


2.9 


0.21 


7.2 


0.630 








4.0 


6.3 


0.46 


8.5 


1.04 







pairs are given in Table II together with the correspond- 
ing values of 5V + . Table II also presents values of y vs and 
i>*; recall: for y$ < y vs , tfi = 0, for y$ — y vs + there 
is a jump of if from zero to if = v*. The most striking 
difference for different (a, a) pairs is in the behavior of 
the Reynolds stress profiles W^(y^) that can be used to 
select the best values of these parameters. 

Clearly, the minimal model with only 4 fit parameters 
cannot fit perfectly the profiles of all the physical quan- 
tities that can be measured. Therefore the actual values 
of a and a should be determined with a choice of the 
characteristics of turbulent boundary layers that we de- 
sire to describe best. Foremost in any modeling should 
be the mean velocity profile which is of crucial impor- 
tance in a wide variety of transport phenomena. Next 
we opt to fit well the profile of the kinetic energy density 
([or, equivalently, the profile of the Reynolds stress tensor 
trace W t (y t )} . Fi gure 7, upper panel, shows the DNS 
profiles of the trace of the Reynold-stress tensor W + {y + ) 
for Re> = 590 (solid lower line) and ReA = 395 (dashed 
lower line). There are no plateau in these plots, meaning 
that these values of ReA are not large enough to have 
a true scale-invariant log-law region. Nevertheless, the 
plots of 

WHy + ) = (l-^r 1 W + (y+) (3.5) 
ReA 

(shown in the same upper panel of Fig. 7) display 
clear plateaus, according to the theoretical prediction for 
ReA — > oo. This means that the decay of W + (y + ) is 
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FIG. 7: Upper panel: DNS profiles of the trace of the 
Reynolds-stress tensor (twice of the total kinetic energy den- 
sities) in two normalization, W + and W*. Solid blue line: 
Re T = 590 and dot-dashed black line - for Re T = 395. 
Lower panel: Comparison of the DNS profile W^(y + ) for 
ReA = 590 (blue solid line) with the results of the sum-MM 
(green dashed line) and the root-MM (red dot-dashed line) 
with constants Eq. (2.13). 



related to the decrease of the momentum flux P(y) and 
that the dimensionless " * " variables, (1.20), that use 
the y-dependent value of the momentum flux P(y) rep- 
resent the asymptotic physics of the wall bounded tur- 
bulent flow, at lower values of ReA than the traditional 
"wall units" (1.18), which are based on the wall value of 
the momentum flux Pq. 

To compare the model prediction with simulational re- 
sults we have to relate W*Ay$) with W^(y + ) in channel 
flows. According to Eqs. (1.18) - (1.20): 
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FIG. 8: Comparison of the DNS Reynolds stress profiles | W£ v \ 
(blue solid lines) for the channel flow with Re^ = 395 and 590 
with the root-MM profiles with constants (2.13)), red dashed 
lines. Upper panel: linear coordinates, lower panel: the same 
plots in the linear-log coordinates. 



and similar Eqs. for the its trace W + (y + ). Figure 7 
shows a peak of W t {y+), W|, ax = W t (y max ) ~ 9.8 at 
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llu 



mal model reproduces the peak in (y* 



18. As one sees from the Tabl. II, the mini- 

*) with an ampli- 
tude of about 8 -j- 8.6 for a < 2. To be specific we choose 
a = 1 in both versions of the minimal model, sum-MM 
and root-MM. With this choice we plot in Fig. 7, lower 
panel, both theoretical profiles, Wj, (y + ) and W* (y + ), 

in comparison with the simulational profile Wo NS (y + ). 
It appears that the root-MM is in better correspondence 
with the simulation than the sum-MM. However, for the 
sake of analytic calculations, the sum-MM is simpler. 
Therefore, again, the choice of the version of MM depends 
on what is more important for a particular application: 
calculational simplicity or accuracy of fit. 



FIG. 9: Comparison of the DNS (solid lines) and root-MM 
results (dashed lines) for partial kinetic energies profiles in 
the channel flow with Re^ = 590. The model parameters are 
taken from (2.13). 



D. Profiles of the Reynolds stress tensor 

In Fig. 8 we present (by solid lines) simulational pro- 
files of the Reynolds stress Wi (y + ) for Re A = 395 and 
He\ = 590 in comparison with the model predictions 
(dashed lines) for the root-MM. The upper panel shows 
the comparison in linear coordinates, the lower panel 
in linear-log coordinates, stressing the buffer layer re- 
gion. In the model profiles we used the values of param- 
eters (2.13), chosen to fit the simulational profiles for the 
mean velocity and the kinetic energy. In other words, in 
comparing the profiles of W xy (y + ) in Fig. 8 no further 
fitting was exercised. Having this in mind, we consider 
the agreement as very encouraging. The only difference 
between the model predictions and the simulational pro- 
files of W+y (y + ) is in a steeper front of the model pro- 
files for y + < 20. This is again because the model does 
not account for the energy transfer that can only flatten 
the front. As already mentioned, even for Re=590 the 
maximum value of the Reynolds stress does not reach it 
asymptotical value |Wi| = 1, as it should in the true 
log-law region. The corresponding comparison for the 
sum-MM looks very similar and is therefore not shown. 

Next, we present in Fig. 9 the simulational and root- 
MM profiles of the diagonal components of the Reynolds- 
stress tensor, W^(y + ) for the channel flow with Re,*, = 
590. Solid lines present simulational profiles, dashed 
lines - the model profiles. The stream-wise and span- 
wise profiles, W£ x (y + ) and W+(y + ), are in good agree- 
ment in the most of the channel, 70 < y + < 470, 
while for the wall- normal component, the model pro- 
file Wy y (y + ) — W^ z (y + ) and differs from the simula- 
tional one. The model also predicts semi-quantitatively 
increase in the streamwise part of the kinetic energy and 
the decrease in the span-wise and wall-normal compo- 
nents in the buffer layer which is observed in Simula- 
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FIG. 10: Profiles of the flatness of the x, y and z compo- 
nents of the turbulent velocity fluctuations, DNS data [15] 
for Re T = 590. Horizontal dotted line shows Gaussian value 
for the flatness, equal to three. 



tions. The physical reason of this is simple: as is well 
known, the energy from the mean flow is transferred only 
to the stream-wise component of the turbulent fluctu- 
ations. Accordingly, in the model one sees the energy 
production term (— 2SW xy ) only in the RHS of equation 
for W xx . The energy redistributes between other compo- 
nents due to "return-to-isotropy" term iy, Eq. (1.14b) 
with the isotropisation frequency oc 1/y. The relative im- 
portance of Iij (in comparison with the energy relaxation 
term) decreases toward the wall due to the viscous contri- 
bution oc 1/y 2 . Accordingly, near the wall only a small 
part of the kinetic energy can be transferred from the 
streamwise to the wall-normal and the span-wise com- 
ponents of the velocity during the relaxation time (that 
oc 1/r). Also, the model describes well the part (about 
50% in the outer layer) of the total kinetic energy that 
contains the streamwise components. 

In the core of the flow (y + > 450) the model gives 
smaller values of all the components , as compared to 
simulations and experiments. This is again because the 
model neglects the energy transfer toward the centcrlinc 
of the channel, where the energy input into turbulence, 
—2SWxy, disappears due to the symmetry reason. 

Also, there is a quantitative disagreement between the 
model and the simulations in the buffer layer. One can 
relate this with the fact that the model neglects the en- 
ergy flux toward the wall, which plays a considerable role 
in the energy balance. The minimal models are local in 
space, but this effect can be effectively accounted for by 
an appropriate choice of the dissipation constants, tak- 
ing a yy > a xx — a zz . We do not propose to take this 
route; in the buffer layer the turbulent flow is strongly 
affected by highly intermittent events (coherent struc- 
tures) connected with the near-wall instabilities of the 
laminal sub-layer. This is confirmed by the very large 
values of the flatness (above 30), as shown in Fig. 10. 



Only for y + > 50 the flatness reaches the Gaussian value 
of 3 and one can successfully utilize various lower-order 
closure model for describing wall bounded flows. 



IV. SUMMARY: STRENGTH AND 
LIMITATIONS OF THE MINIMAL MODEL 

The minimal model as formulated in this paper is a 
version of the algebraic Reynolds stress models. Its aim 
is to describe, for wall bounded turbulent flows, the 
profile of mean flow and the statistics of turbulence on 
the level of simultaneous, one-point, second-order veloc- 
ity correlation functions, i.e. the entries of the Reynolds- 
stress tensor Wij. The model was developed explicitly 
for plain geometry, including a wide variety of turbulent 
flows, like channel and plain Couette flows, to some ex- 
tent fluid flows over inclined planes under gravity (mod- 
elling river flows), atmospheric turbulent boundary layers 
over flat planes and, in the limit of large Reynolds num- 
bers, many other turbulent flows, including pipe, circular 
Couette flows, etc. 

In developing a simple model one needs to decide what 
are the physically important aspects of the flow statistics, 
those which determine the mean-flow and the turbulent 
transport phenomena. The choice of the Reynolds-stress 
approach was dictated by the decision to emphasize the 
accurate description of V(y)- and Wij(y)- profiles. The 
main criteria in constructing the model were simplicity, 
physical transparency, and maximal analytical tractabil- 
ity of the resulting model. That is why we took lib- 
erty to ignore the spatial energy flux, and, thanks to 
the plain geometry, to estimate the spatial derivatives 
and the outer scale of turbulence using the distance to 
the wall y. The same motivations led to choosing the 
simplest linear Rotta approximation of the "Return to 
isotropy term" [19] and the simplest dimensional form of 
the nonlinear term for energy flux down the scales, also 
in agreement with [19]. 

By proper parametrization the number of fit parame- 
ters was reduced from twelve to four. Two of these, a 
and a are responsible for the viscous dissipation of the 
diagonal, Wu , and the off-diagonal, W xy , components 
of the Reynolds-stress tensor. The other two parameters 
- b and b - control the nonlinear relaxation of Wu and 
W X y It appears that one cannot decrease the number 
of fit parameters further with impunity. The outer layer 
parameters b = 0.256 and b — 0.500 where chosen to 
describe the observed constant values of von-Karman in 
the log-law (2.10) and the asymptotic level of the den- 
sity of kinetic energy. The viscous layer parameters were 
chosen to describe the observed values of the intersection 
C in the von-Karman log-law (2.10) and the peak of the 
kinetic energy in the buffer sub-layer. The resulting set 
of 5 equations for the mean shear S(y), Reynolds stress 
W xy and W xx , W yy , W zz with just four fit parameters is 
referred to as the minimal model. 

As demonstrated in Sec. Ill the minimal model with 



the given set (2.13) of four parameters describes five func- 
tions: 

• the mean velocity profile V(y) is describe with ac- 
curacy of ~ 1% - almost throughout the channel 
(except of small velocity defect in the core of the 
flow), cf. Fig. 6; 

• the Reynolds stress profile W xy (y) is also described 
with accuracy of few percents (except in the viscous 
layer y + < 5 in which W xy does not contribute to 
the mechanical balance), cf. Fig. 8; 

• the total kinetic energy profile \ W(y) is reproduced 
with reasonable (semi-quantitative) accuracy, in- 
cluding the position and width of its peak in the 
buffer sub-layer, cf. Fig. 7; 

• The profiles of the partial kinetic energies, 
^W xx (y), \W yv {y) and \W zz {y), are reproduced, 
see Fig. 9, including the simple \-\-\ distribu- 
tion in the asymptotic outer region. This distribu- 
tion is supported by recent experimental, DNS and 
LES data, as shown in Figs. 3, 4, and 5. 

We consider all this as a good support of the minimal 
model; too much data is being reproduced to be an ac- 
cident. It appears that the minimal model takes into ac- 
count the essential physics almost throughout the chan- 
nel flow. 

On the other hand, one should accept that such a sim- 
ple model cannot pretend to describe all the aspects of 
the turbulent statistics in wall bounded flows. For exam- 
ple, the minimal model ignores the quasi-two dimensional 
character of turbulence and the existence of coherent 
structures in the very vicinity of the wall. The minimal 
model does not attempt to take into account many-point 
and high-order turbulent statistics, including three-point 
velocity correlation functions and pressure-velocity cor- 
relations, responsible for the spatial energy flux and for 
the isotropization of turbulence. Finally, our choice of 
dissipation term definitely contradicts to the near-wall 
expansion, (and see Sec. 11.7.4 of [5]), in disagreement 
with various known improvements of [19]. We propose 
that all this is a reasonable price for the simplicity and 
transparency of the minimal model, which is constructed 
with emphasis on the fundamental characteristics V(y) 
and Wij(y) which are crucial for most applications. 

We trust that a proper generalization of the minimal 
model will be found useful in the futures in studies of 
more complicated turbulent flows, laden with heavy par- 
ticles, bubbles, etc. 
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APPENDIX A: VALIDATION OF THE 
ITERATIVE PROCEDURE 

To see, how the iterative procedure described in 
Sec. II B 2 works, we plotted in Fig. 11 iterative profiles 
of the turbulent velocity v^(y$) for n = 1, ... 5 together 
with the (numerical) solutions of Eq. (2.17a), v+(y$) (the 
thick solid line) and v_(y^) (dot-dashed curve). The hor- 
izontal straight line presents the asymptotic value v^. 
The critical point {vl, yl s } is shown by a black circle. 
Our analysis shows (and see also Fig. 11), that already 
the simple Eq. (2.22a) gives the relative accuracy (with 
respect to v^) better than 1% for y$ > 30. The second it- 
eration works with this accuracy in wider region y$ > 10, 
the third iteration gives 1% accuracy for y$ ss 5, which is 
about the critical value y\ s ss 4.8. Unexpectedly, the ap- 
proximate solutions work even below the y* s) where exact 
solution is v$ — 0. One observes with increasing n the 
widening of the region, in which practically indistin- 
guishable from zero. The overall conclusion from these 
observations is that already the fist few iterations give a 
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very good accuracy for all practical purposes, and very often one can use only the first or the second iteration. 
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